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ABSTRACT 

We derive similarity solutions for the expansion of negative initial density pertur- 
bations 6M/M cx (s > 0) in an Einstein-de Sitter universe filled with collisional 
baryonic gas and coUisionless matter. The solutions are obtained for planar, cylindrical, 
and spherical perturbations. For steep perturbations, central cavities surrounded by 
over-dense shells expanding as t2(«+i)/(3s) develop in both matter components. We also 
consider the case when baryonic shells are driven by internal pressure. Before redshift 
~ 8 — 10 baryonic shells cool by inverse Compton scattering with the cosmic back- 
ground radiation and can form bound objects of mass ~ lO^^Q^^^"^ Mq{1 + z)^'^/'^ . Hot 



shocks of radius R cause a Compton y-distortion of order y 
where fp is the volume filling factor of these shells. 



10-5r2b(i?/10Mpc)2/p, 
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1 INTRODUCTION 

Gravitational collapse of matter in initially over-dense re- 
gions leads to the formation of dense bound objects in the 
universe. Star formation is efficient in these objects and so 
we expect them to harbor most of the luminous matter in 
the universe. In the linear regime, gravitational instability 
simply amplifies the density fiuctuations by a universal time 
dependent factor. Thus, an initially under-dense region re- 
mains under-dense and cannot lead to the formation of dense 
structures. In the nonlinear regime, however, the evolution 
of under-dense regions depends strongly on the profile of 
its density perturbation. In an isolated under-dense region, 
the density remains below the background if the profile is 
shallow. If the profile is steep then expanding inner shells 
can catch up with shells farther out causing the appear- 
ance of caustics in the coUisionless matter component and 
shocks in the baryonic gas (e.g., Fillmore & Goldreich 1984, 
Bertschinger 1985). Caustics and shocks are very dense and 
can fragment under their self-gravity to form bound struc- 
tures. 

In a generic density field, a large scale under- or over- 
dense region contains small scale substructure of density 
fluctuations. In a large scale over-density, substructures 
merge under the large scale gravitational field and become 
part of a large object typically of the mass of the large 
scale region. Dynamics of under-dense regions can be more 
complex, depending on the details of the large scale den- 
sity profile. Over-densities develop into virialized bound ob- 
jects which are easily identified in observations and N-body 



simulations, by means of the friends-of-friends algorithm for 
example. The fate of under-dense regions however is less cer- 
tain and strongly depends on the initial conditions. Because 
of the complexity in studying the dynamics of matter in 
under-dense regions, analytic solutions for special symmet- 
ric negative perturbations are particularly important. Most 
authors focused on symmetric initial perturbation of the 
form 5M/M oc r"". In a matter dominated flat universe the 
gravitational evolution of this proflle leads to self-similarity, 
which considerably simplifies the problem. Another simpli- 
fying assumption, that is often made for evolution of under- 
dense system, is that of spherical symmetry. To see the phys- 
ical motivation for this assumption consider an ellipsoidal 
top-hat perturbation. Since the perturbing force (accelera- 
tion for negative and deceleration for positive perturbation) 
is stronger along the minor axis, the asymmetry of a void 
will decrease. Ryden (1994) has investigated the evolution of 
under-dense coUisionless axisymmetric systems with the ini- 
tial perturbation of the form SM/M oc r~''f{0), where O is 
an angle with the symmetry axis, and found that the emerg- 
ing void is nearly spherical. Fillmore & Goldreich (1984) 
have found similarity solutions for spherical perturbations in 
coUisionless matter with s < 3, while Bertschinger(1985) ob- 
tained solutions for s = 3 (uncompensated hole) for both gas 
and the coUisionless matter. However, under certain condi- 
tions, the system becomes unstable and spherical symmetry 
might be destroyed. Several authors have developed different 
approximations to determine the conditions for instability 
and to evaluate the sizes of the formed objects (Bertschinger 
1983; Vishniac 1983; Hwang, Vishniac and Shapiro 1989; 
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White and Ostriker 1990). In the special case of a compen- 
sated void (i.e. void surrounded by an over-dense shell whose 
mass is equal to the mass deficit in the void) , the dense shell 
typically fragments into large massive clumps (White and 
Ostriker 1990). In this paper wc investigate the self-similar 
evolution of negative perturbations for all positive values of 
s. Our main focus is spherically symmetric perturbations, 
but, in order to see how the evolution depends on geometry, 
we also present the main features of planar and cylindrical 
perturbations. 

The paper is outlined as follows. In §2 wc write the 
equations of motion and transform them into the appropri- 
ate self-similar form. In §3 we derive the asymptotic behavior 
of the equations for various values of s, and the adiabatic in- 
dexes, 7. In §4 the qualitative features of the numerical solu- 
tion arc described. In §5 we check the possibility of structure 
formation in coUisionless and baryonic shells. In §6 we con- 
sider evolution of shells, whose expansion is dominated by 
an internal pressure source. In §7 we present the summary 
of our results. 



2 THE EQUATIONS 

Wo write the Newtonian equations of motion governing the 
adiabatic evolution of symmetric perturbations in a coUi- 
sional fluid (gas) of adiabatic index 7. We assume that 
the expansion scale factor of the universe is a{t) oc t^^"^, 
where t is the age of the Universe, the Hubble function 
is H{t) — 2/{3t), the total background density is Pc = 
3H'^/{8TrG) = l/{6nGt^), and the ratio of the mean col- 
lisional baryonic density to the critical density pc is denoted 
by Qb- Although self-similar evolution exists for any flh < 1, 
in this paper we restrict the analysis to either fib = 1, or 
fib < 1- 

Denote by r and v = dr/dt the physical position 
and velocity of a gas shell, where r = is the symmetry 
center of the perturbation. Further, let p{r,t) and p{r,t) 
be the gas density and pressure at r. Wc write the mass 
within a distance r from the symmetry center as m{r, t) = 
x^''^^ p{x,t)dx, where n = 1,2, and 3 refer, respectively, 
to planar, cylindrical, and spherical perturbations. The mass 
within a flxed shell varies with time like m ~ j-2(3-n)/3^ 
cause of the Hubble expansion along 3 — n of the axes. At 
this stage we neglect such physical effects as viscosity, ther- 
mal conduction, cooling etc., so the equations of motion of 
coUisional fluid are, the continuity equation. 



2(3-n) 
pt 3 



2(3-r.) 



-t^^pr'-"a^(r"-'v) 



At 



the Euler equation, 

dv 2 3 — n r _ drP AirGmx 
di ~ 9 n ¥ ~ p r"-i 

the adiabatic condition, 
^(pp-^) = 0. 



(1) 



(2) 



(3) 



and the relation. 



(-\ n — 1 

^T-m = r p . 



(4) 



In equation (2), rux stands for total (coUisional and col- 

lisionless) mass. 

There are several similar techniques, which we do not 
describe here in detail, for calculating the coUisionless mass 
profile (Bertschinger 1985; Fillmore & Goldreich 1984). The 
main idea is that, since the motion is self-similar, knowing 
the orbit of a single particle one can calculate the mass pro- 
file and vice versa. Given the initial guess for mass profile, 
after iterative calculation of mass profile and orbit one ar- 
rives to the self-similar solution. 

The initial conditions leading to self-similar expansion 
are specified at an early time close to zero, U, as 

6M 

H 

2 

% 
0, 



,r » ro 



(5) 



v{r,ti) = 
P{r,ti) = 



(6) 
(7) 



where s > and 6M/M is the mean density contrast in- 
terior to r. The Einstein-De Sitter Universe and the initial 
conditions are completely scale-free. The only characteristic 
length in the evolution of the perturbation is the scale of 
non-linearity r» (t) . For self-similar collapse of positive per- 
turbations the turn-around radius is usually used as the scale 
of non-linearity. For a negative perturbation, r* can be de- 
fined at any time as the radius interior to which the mean 
density contrast has a certain fixed value. This means that 



r* oc r 



2s 



3(s + l) 



(8) 



This is also the way the turnaround radius in positive pertur- 
bation depends on time. The choice of proportionally factor 

in (8) is arbitrary. Here it is chosen such that a particle with 
initial radius reach r, at {65M/M)^'^^ ^ti, where the fac- 
tor of 6 is arbitrary. For this choice of r, the mean density 
contrast interior to r* is approximately —0.09 in all symme- 
tries. 

A consequence of self-similarity is that the partial differ- 
ential equations of motion can be transformed into ordinary 
differential equations. This is done by working with A = r/r* 
and the dimensionless ffuid variables. 



v{r,t) = ^y(A) 

p(r,t) = QbPcDiX) 

p{r,t) = fibPc(j)'i^(A) 

m{r,t) = ifibPcr"M(A) . 



(9) 
(10) 

(11) 

(12) 



Expressed in terms of these variables, the equations (1- 
4) become, respectively. 



(V -aX)D' + + V' -f) D = , 



(13) 
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(^^^-^yV-a\) = 2{a-2 + ^) , (15) 

M' = 3A""^£) , (16) 

where the prime symbol denotes derivatives with respect to 
A. 

Self-similarity implies that the shock appears (if it does) 
at fixed A = As = rs/r*, so the physical radius of the shock 
Ts oc t" and its non-dimensional speed is (r*/t)~^(drs/dt) = 
a As. At the surface of the shock the fluid variables satisfy 
the jump conditions obtained from mass, momentum, and 
energy conservation (cf. Spitzer 1978, §10.2): 



the energy inside the shock must increase as it sweeps more 
particles. On the other hand, self-similarity implies that the 
total energy within As varies with time as which is a 

decreasing function of time for a < 4/5, i.e., s < 5. To sat- 
isfy both conditions the energy within As must bo negative. 
The only way to construct such system is by placing a bound 
object (i.e., negative energy) at the center. However, since 
there is no way to make this object grow with time, the evo- 
lution of the perturbation never becomes self-similar. The 
limiting case of s 5, which corresponds to a shock in a 
homogeneous medium (compensated void), was investigated 
by Bertschinger (1983). The corresponding upper limits for 
cyhndrical and planar perturbations are s fs 3.5 and s = 2, 
respectively. 



2.2 Boundary conditions 



f 

v'_ 


P- 7-1 
P+ 7 + 1 ' 


(17) 


P+ + P+V+ 




(18) 


V- 




(19) 


v+ 


; , drs 


(20) 



where the superscripts of the minus and plus signs refer to 
pre- and post-shock quantities and v' is the velocity rela- 
tively to the shock position. 

In terms of the non-dimensional fluid variables we ob- 
tain 

V+ = aXs + ^^{V- -aXs) , (21) 
7 + 1 

D+ = 1±1d- , (22) 
7-1 

P+ = -^D-{V- -aXsf . (23) 

7 -I- 1 

Outside the shock, at A > As, the pressure vanishes and the 
prc-shock fluid variables can be found by solving the equa- 
tions (13-16) with zero pressure. Analytic solutions for zero 
pressure exist for planar and spherical geometries, but not 
for cylindrical (Zel'dovich 1970, Peebles 1980, Bertschinger 
1985). 



2.1 The existence of shocks and the value of s 

Lower and upper limits on s exist for a shocked expansion of 
an isolated negative cosmological perturbation. A necessary 
condition for the formation of a shock is that the perturba- 
tion is steep enough so that inner shells expand faster than 
outer shells. An inspection of the equations of motion shows 
that this happens only if (a — l)a < f that is s > 2 
for spherical, s > 1.54 for cylindrical, and s > 1 for planar 
symmetry. This lower limit disappears if an external energy 
source is placed at the center of the perturbation. 

An upper limit on s also exists. To see this, consider 
the case of a spherical perturbation. If the perturbation is 
negative, the total energy of a particle is positive and so 



Solving the self-similar equations (13-16) requires boundary 
conditions on the fluid variables. We will see in the next sec- 
tion that shocked expansion is associated with a dense shell 
whose outer boundary is the shock. Outside the shock the 
pressure is zero and the equations can be solved by standard 
numerical techniques. The jump conditions (21-23) provide 
post-shock fluid variables. Given the post-shock fluid vari- 
ables, the equations (13-16) can be integrated inward to 
obtain the fluid variables inside the shock. 

So far it seems that one can introduce the shock at any 
arbitrary radius. However, inward integration of the equa- 
tions from an assumed shock position always end at a sin- 
gular point, Ao, where at least one of the fluid variables 
becomes inflnite. Since there is no natural way of eliminat- 
ing this singularity, this must be the inner boundary of the 
baryonic shell, Mgas (Ao) = 0. If the total energy of the sys- 
tem is to be conserved, the pressure at the inner boundary 
of the system must be zero, -P(Ao) = 0, otherwise, energy 
will continuously be injected into the system. It turns out 
that there is at most one value of As , for which the solution 
satisfies -P(Ao) = 0. Smaller values of As give positive mass 
and pressure at Ao, while larger values give zero mass and 
positive pressure. 



3 ASYMPTOTIC BEHAVIOR NEAR THE 
INNER BOUNDARY OF THE SHELL 

3.1 Baryonic matter 

To obtain an asymptotic behavior for the collisional matter 
we expand the fluid variables near Ao, as 



M(A) 


= Mo{X-Xor 


D{X) 


= Do{X-Xof 


P{X) 


= PoiX-Xor 


V{X) 


= aXo + Vo{X- 



XoY 



By substituting these expressions in the equations of 
motion we flnd that f = 1 in all cases(cf. Ostriker & McKee 
1988), while 6, and rj depend on s and 7. The asymptotic 
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exponents of the fluid variables in the mixed and purely col- 
lisional matter are the same. This differs from the collapse 
case where the collisional component can change the asymp- 
totic exponents of the collisional fluid variables near the cen- 
ter (Chuzhoy & Nusser 2000). Table 1 lists the asymptotic 
exponents for the following three cases: 

/. Expansion without a shock, Ao — 0: As we have men- 
tioned in the previous section, there are no solutions with a 
shock for (a - l)a > with P(Ao) = 0. The expansion 

then develops without a shock and the pressure vanishes 
everywhere. Substituting the asymptotic form for the fluid 
variable with Ao = 0, in the equation of motion yields the 
asymptotic exponents in this case as listed in Table 1. Near 
A = 0, the mass vanishes and the velocity grows linearly as 
VqA, where Vo is independent of s. 

IT. Shocked expansion with P(Ao) = 0; This occurs if 
(q - l)a < §2^, i.e., if s > 2 for spherical, s > 1.54 
for cylindrical and s > 1. For these values of s the exponent 
5 as given in table 1 is always negative and so the density 
diverges at Aq. 

///. Shocked expansion with J-'(Ao) > 0; This occurs when 
energy can be injected from the center of the perturbation 
(e.g., SN energy injection) allowing -P(Ao) — const > at 
the inner shell boundary. The asymptotic behavior for this 
case is also listed in the table. 



3.2 Collisionless matter 

After going through first shell-crossing, a collisionless parti- 
cle starts its converging oscillatory motion around some A^, 
inside the dense shell (Fillmore & Goldreich 1984). At each 
turnaround point there is a caustic peak, where the density 
diverges. To obtain the asymptotic solution at the caustic 
we again expand the variables around the peak location 

M(A) = Mo(A-Ao)^ 
D{X) = Do{X-Xof 
P(A) = 

y(A) = aAo + VoiX - XoY . 

Substituting into the continuity, Euler and mass equations, 
we obtain /i = 1/2, 5 = -1/2, ly = 1/2 and Vq = 2(1 - 
Q)aAo. Negative and positive Vo correspond to the incoming 
and out-coming streams at the caustic. Note that, unlike the 
collisional case, the asymptotic exponents depend neither on 
symmetry nor value of s. 



4 NUMERICAL SOLUTIONS 

In order to obtain the detailed structure of the baryonic shell 
we numerically integrate the self-similar equations of motion 
(13-16). The location of the shock, As, is sought by demand- 
ing vanishing pressure at the inner boundary of the bary- 
onic shell. Integrating the equations inwards from arbitrary 
As always leads to a singular point, where the derivatives of 
the fluid variables diverge. If the assumed As is larger than 
the desired value, then at this point pressure is positive and 



Table 1. Asymptotic constants, Vb, 5, ?? collisional fluid variables 



in a flat 


universe with no 


collisionless matter. 






Ao =0 


Ao > 


Ao > 




NO SHOCK 


P(Ao) = 


P(Ao) > 


S 




n-2-3(37-4) 


s(37-4)-H2 




i\/n+'24-{s+4)y/n 


2-Hn(7-l)-H3(37-4) 


2-Hn7-Hs(37-4) 




5-1-3 


5-1-1 


S + 1 


V 


P = 


S+1 








2(n-2-|-7(l-n)-|-2s(7-2)) 




Vo 




37s 




mass is 


zero, while if it 


is smaller then the mass 


is also pos- 



itive. Thus it is possible to find As iteratively starting from 
the initial guess and varying it in accordance with the re- 
sults of the integration. We will show numerical results for 
only two cases fib = 1 and f2b ^ 1- To find the solution for 
fib <SC 1 we need the mass profile for the purely collision- 
less collapse. This is found using the method of Fillmore&: 
Goldreich(1984) and Bertschinger(1985). 

Figures 1-2 show the shock location and thickness as a 
function of s. Figure 3 shows the fiuid variables as a func- 
tion of A for several values of s. All the plots were made 
for a spherical expansion with 7 = 5/3. All the numerical 
solutions we found conform near the inner boundary of the 
shell to the results of asymptotic analysis. The location of 
the shock (fig. 1) depends strongly on s but is nearly in- 
dependent of both of 7 and fib. The reason for this weak 
dependence on 7 is that the thermal energy is small to the 
kinetic and gravitational energies. The thermal energy and 
hence the width of the shell increase with s (fig. 2), from zero 
to a few per cent of the void radius. The baryonic shell lies 
in the collisionless shell. Because the collisionless particles 
are not slowed down as they enter the shell, the collisionless 
shell is typically several times wider than the baryonic shell. 



5 FRAGMENTATION 
5.1 Baryonic matter 

Wo address the question whether the shells can fragment 

into bound objects, such as galaxies. 

For the expansion of an object to halt, its gravitational 
energy must be larger than the sum of the thermal and the 
kinetic energy by some finite factor, i.e. 



Eg + Cf{Et, + Et) <0 (24) 



© 0000 RAS, MNRAS 000, 000-000 



self-similar expansion 5 




2.5 3 



3.5 
s 



4.5 5 



Figure 1. The shock location As as a function of s. The dashed 
and the solid line correspond, respectively to Ob = 1 and fit *S 1- 




4.5 5 



Figure 2. Shell width AA/As as a function of s. The solid and 
the dashed line correspond, respectively to Q), = 1 and *S 1- 



To check whether an object satisfying the above equa- 
tion can be found, we calculate the energy of a disk placed 
on the interior boundary of the shell (where the density is 
highest and the temperature is lowest). 

Let the disk be of radius Z = r,L and height h = r*H. 



The kinetic energy due to the shear velocity is found to 



be 



(25) 



2 

0.75 

0.65 

0.55 

0.45 
10-^ 




10"- 



H \ 1 \ \ 1 h 



. ...U- 




■=*p-| ^ V- 



10^ 

0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 

"o 

Figure 3. The fluid variables (density, velocity, pressure and ther- 
mal energy) for the case of spherical symmetry and 7 = 5/3. The 
solid and dotted lines correspond to ilt = 1 with s = 3 and s = 4 
respectively , dashed and dash-dotted lines correspond to f2b ^ 1 
with s = 3 and s = 4. 



3, respectively and ma « f^Mnbpcf"* is the mass of the 
disk. 

The gravitational energy is obtained using the thin shell 
approximation * - 



The thermal energy ■ 



2,3M2(A„+i/) 



\4 



(26) 



where C„ = 4/9, (2/9 + a^/2), and forn = 1, 2, and 



* For all the disks satisfying (24), that we found, H <^ L. Thus, 
since the correction for finite thickness of the disk lowers the grav- 
itational energy, the assumption of a thin disk is justified. 
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'Ao 

The kinetic energy due to the radial velocity ■ 



(27) 



-Bk2 = / irfpiy- < V >Y/2dr ■■ 



D{V- < V >y/2d\ 



(28) 



(29) 



where < V > is an average radial velocity of the disk. Prom 
the asymptotic analysis follows that £k2/-E't oc H, so for 
J? <§; 1 we can neglect E\^2- 

For large values of L the kinetic energy term, E]^ oc L*, 
would be dominant, while for small L the thermal energy, 
Et oc L^, is the only significant term. Thus the disk can 
be unstable only for intermediate values of L, where the 
gravitational term, Eg oc L^, is dominant. From eq. (24) we 
can factor an equation of second degree in L, whose roots 
determine the limits of instability for a given thickness of 
the disk, H. The condition for the existence of the roots is 



i6ab(7-i)-rr > 3Vac| / 



Xo+H 



PdX 



(30) 



5.1.1 Fragmentation without cooling 

The asymptotic solution implies that the term on the l.h.s. 
of eq. 30 grows as H^C^+i) ^^le right as and thus 

in the limit — » the inequality is satisfied if 5 < — i . 

For 7 = 5/3, 5 is less than — i for all planar shocks 
and for cylindrical shocks with s > 10/3, which are conse- 
quently gravitationally unstable even for <^ 1. For spher- 
ical shocks 5 is less than — i only in case the shell pushed by 
an internal pressure source with expansion rate a < 8/9. For 
spherical shocks with P{Xo) = the numerical integration 
shows that for f2b = 1 and C/ = 1 eq. (30) is satisfied only 
for s >~ 3.5 and not at all if Qb < 0.8. 



5.1.2 Fragmentation with cooling 

From the numerical calculation we find that for a disk on a 
spherical shell with zero internal pressure to satisfy eq.(30) 
its thermal energy must fall at least by factor Q,^^ . To check 
when this is indeed possible we must introduce a physical 
scale into the problem. 



IS 

T: 



The physical temperature of the shock for ionized gas 
Pa 0.6m„ 



DaXi kB 



'■p_ ( ' 

V t ) 



(31) 



For spherical symmetry the value of Ps/DgAg can be 
approximated by (s — 2)/ 10^ for shell driven by negative 
density perturbation (fig. 4) or by 3(a — 2/3)^/16 for shell 
driven by internal pressure. Thus 



T«2(s-2)loVio(l + 2) 



or 



{s-2)/a 



K 



T^4{a- 2/3)^loV?o(l + zf-'^'^^K, 



(32) 



(33) 



where rso = lOno Mpc is the present radius of the 
shock. 

At temperatures higher than ~ lO^K and densities 

not much higher than the background values the inverse 
Compton scattering is the dominant cooling process. At red- 
shifts > 8 the cooling time tcooi = 1-6 • 10^^(1 + z)~*Yt 
(Bertschinger 1983) is less than the age of the universe. This 
means that the gas shock cools quickly after the shock and 
the average temperature of the shell remains at constant 
value of order of 10^ ii". Thus the approximate condition for 
fragmentation is for the shock temperature to be higher than 

When the collapse depends on cooling the mass of the 
formed objects is of order of Jeans mass at T ~ W^K, i.e., 

M ~ lO"n^^''^(l+2:)"^''^M0. Since the mass accretion rate 
on the shell is proportional to (1 -|- z)'"''^^"'^'^, the number 
of objects of the mass M formed by the shell is proportional 

to N{M) oc M'^"-''\ 

It should be noted that different equations of motion 
(Appendix A) are needed to describe the shell that loses its 
thermal energy during fraction of its dynamic time (snow- 
plow). The main effect of the cooling is to move the shock 
inwards thus decreasing the swept mass and the shock tem- 
perature. However, since in this problem the thermal energy 
constitutes only a small fraction of the total, these differen- 
cies are generally insignificant. 



5.2 Collisionless matter 

We now address the stability condition of the collision- 
less shell. The kinetic and gravitational energy terms are 
the same as in the baryonic case. There axe two stream 
fiows near Ao (cf, section 3.2), so the "thermal" energy 
resulting from the velocity dispersion in the flow is E^ = 



Svrf 

case IS 



f /x„ D{V - aXfdX. The analog of eq. 30 for this 



32^ > 3V^CyCn / ° D{V - aXfdX 

J An 



(34) 



Near the inner boundary both the left and the right term 
grow as H^^^ and thus the fact whether the last inequality 
is satisfied depends on the choice oi Cf. However, using just 
the lower constraint Cf > 1, we found that for small H 
it is never satisfied in the case of spherical symmetry. For 
large H the energy analyses must give the same results for 
pure coUisional {Qh = 1 and Qc = 0) and pure collisionless 
(Jlfc — and Oc — 1) models. Thus we expect the bound 
structures to form for all s in planar shells, for s >^ 3 in 
cylindrical shells, and possibly for s >~ 3.5 in spherical 
shells. 
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3.5 




2 2.5 3 3.5 4 4.5 5 
s 

Figure 4. Non-dimensional shock temperature as a function of s 
(eq. 31). The solid and the dashed line correspond, respectively, 
to Hb = 1 e'lid Qb 1- 



6 EXPLOSIVE AMPLIFICATION 

The explosive amplification model suggesting that galaxies 
can originate from scries of explosions was proposed by Os- 
triker& Cowie(1981) and later developed by other authors. 
Vishniac, Ostriker & Bertschinger (1985) showed that an 
energy injection from supernova explosions in the range 
10^^ — 10®^ could produce a void reaching ~ lOMpc. The 
cooling of the shell surrounding the void would then load to 
the formation of galaxies. However, since the cooling pro- 
ceeds chiefiy through inverse Compton scattering, it was 
argued that a large number of such explosions would be 
incompatible with the observational limits on the Comp- 
ton y-distortion parameter (Levin, Freese & Spergel 1992). 
Later Miranda&; Opher(1997) showed that if the energy was 
injected in several rather than in a single cycle, then the dis- 
tortion would be several times smaller. In their calculation 
the current observational limit y <^1.5 ■ 10^'' (Fixen et al. 
1996) is compatible with voids of present radii 7? ~ 10 — 30 
Mpc, the filling factor fp ~ 0.07-0.30, and fib = 0.1. In this 
section we explain how the distortion produced by shocks of 
similar size might be still lower by one or two orders of mag- 
nitude. 

In the previous section we have shown that before 
z ^ 8 — 10 large baryonic shells become gravitationally un- 
stable. The collapse of one object creates a pressure gradient 
causing the surrounding gas to expand and thus creates a 
feedback against further collapse. Therefore we expect that 
fragmentation proceeds continuously rather than in cycles. 
Depending on the star formation efficiency in the collapsed 
objects and the fraction of stars turned to supernovae (SN), 
the dynamics of the shell might become dominated by the 
energy input from the SN explosions. By using the approxi- 
mation Linj oc msf , where Linj is the energy injection rate 
and nie is the mass accretion rate on the shell, we get the 
expansion rate of the shell rs oc t", where a = 1 + q/2. Note 



that if mass to energy conversion efficiency doesn't change 
drastically over time, then q m and awl. 

If the shell expansion is dominated by internal pressure 
then the collapsing objects, whose surface area is decreas- 
ing, would be eventually overtaken by the shell and, unless 
Jib ~ 1, remain inside. Since the interior of the shell is empty 
of gas, only a small fraction of the SN energy is turned into 
a thermal energy, making the CBR distortion in this sce- 
nario much smaller than in previous models (see Appendix 
B). Similarly, because less energy is lost through cooling, 
larger voids can be produced for the same mass to energy 
conversion efficiency. In this case the radius of the bary- 
onic shefi is w 10/1-^(^5^)^/^(1 + z)'^''^ Mpc, where 
e = Einj /mshM(? ■ 



7 SUMMARY 

For the gaussian initial density field with scale free power 
spectrum the rms density perturbation is < {^^Y 

where — 3 < n < 4. It is therefore reasonable to 
identify s = {n + 3)/2, from which follows < s < 3.5. 
Thus in case n > 1 we might expect that the universe con- 
tains a large number of spherical shells. We have shown in 
the paper that the evolution of spherical perturbations with 
s > 2 results in formation of overdense regions, which, unless 
(s — 2) -C 1 may contain a large fraction of matter. 

A distinctive property of structure formation in the shell 
is the segregation between the baryons and the coUisionlcss 
matter. The energy analysis shows that for s < 3.5 (i.e. 
n < 4), which follows from the gaussian density field, the 
coUisionless matter cannot form bound objects, while at high 
redshifts cooling allows formation of the baryonic objects on 
the appropriate scales. Similarly, the explosive model (§6) 
predicts that the baryonic shell is driven by pressure ahead 
of the coUisionless and fragments in an independent way. 
This is entirely different from the evolution of positive per- 
turbations, where gas must fall into the potential well formed 
by coUisionless matter. 

The energy analysis (§5.1.2) shows that the mass of the 
bound objects formed from the baryonic shell is of order of 
Jeans mass at T ~ IG^K - M ~ 10"n-^/^(l -|- z)-^''^Mq. 
This is again different from the evolution of positive pertur- 
bation, which predicts the Jeans mass to be also dependent 
on the spectrum of perturbations. 

An observational limit on the scale and rmmbor of the 
shells can be obtained from the amplitude of Compton y- 
distortion. Both for the shell driven by the negative density 
perturbation and for the shell driven by internal pressure the 
predicted amplitude is of order 10-^^(-^)^/f, where /f 
is the present filling factor (see Appendix B for the details 
of the calculation). 

The current upper limit y = 1.5- 10"* is thus consistent 
with shells, whose present radius is below ~ 30 Mpc and the 
filling factor is of order of unity. 
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APPENDIX A: EXPANSION OF A THIN MOMENTUM CONSERVING SHELL (SNOW PLOW) 

If the cooUng time is much shorter than the dynamical time of the system the shocked gas would be compressed into a thin 
shell, whose expansion can still be described in self-similar form (r^ = raot"). Given the initial density profile and energy 
injection rate, the constants a and rso can be determined using momentum conservation, i.e., 

where fp is the force of the internal pressure, p and v are the density and the velocity of the gas just outside the shock, and 
m and ttIc are, respectively, the mgisses of baryonic and coUisionless components. In the dimensionless form (Al) is expressed 
as, 

2(3 -n) 27W(aA - l/)A"-2 K, , . 

6a(a - 1 + n(a - 2/3)) = - C- + L__^ + ^, (A2) 

where C = 2 for fib ^ 1 (if the baryonic shell is ahead of the coUisionless) and C = 1 for fib = 1. 

If the expansion is caused by instantaneous energy injection then a is found by substituting M = 3A/n, D = 1, V = 2A/3 
and = in (A2). The result is 



3 + 4,; + ^{21- Tm - i2C[u + l))/n 

a = r . (A3) 

6(n+l) ^ 

For rib = 1 a is equal respectively to 0.797 for spherical symmetry, 0.853 for cylindrical and 1 for planar symmetry. For 

fib ^ 1 and the baryonic shell is ahead of the coUisionless a « 2/3 + fib/2n. 

If on the other hand the expansion is caused by negative density perturbation or a continuous energy injection then while 
a is the same as in the adiabatic case where cooling is negligible. This is because the thermal energy is only a small fraction 
of the total energy anyways. Cooling however causes a slight decrease in Tso- 



APPENDIX B: SHOCK IMPRINT ON THE COSMIC BACKGROUND RADIATION 

The spectral distortion of the cosmic background radiation (GBR) is given by 

where and Te are the electron density and temperature and ctt is the Thompson cross section for photon electron scattering. 

Averaging over a spherical volume gives < dy/dl >= arEth/ {'^nmeCprt), where Eth is the thermal energy of the entire 
sheU. 

Prior to z ^8 the shell cools quickly so its thermal energy may be approximated by 

Jo ' ' ^o..dt=—-kBT.t^ooi = {— (B2) 

where p is the mean mass per particle. At later redsliifts the thermal energy is Eth ~ J^" Aixpr^dr ~ ( '^ct'^ )(^^'^ jx' P'^^)- 
The filling factor of the shells grows as f(t) = fpt^"'^, where fp is the final filling factor, thus the final distortion is. 

Both for spherical shells driven by the internal pressure (4/5 < a < 2) and shells driven by the negative density 
perturbation y ^lO^^/Friofib, where fp and lOrio Mpc are the present filling factor and radius of the shells . 
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